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Novel Phase Transition Near the Quantum Critical Point in the Filled-Skutterudite 

Compound CeOs 4 Sbi 2 : an Sb-NQR Study 
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We report on a novel phase transition at T — 0.9 K in the Ce-based filled-skutterudite 
compound CeOs4Sbi2 via measurements of the nuclear-spin lattice relaxation rate 1/Ti and 
nuclear quadrupole resonance (NQR) spectrum of Sb nuclei. The temperature (T) dependence 
of 1/Ti behaves as if approaching closely an antiferromagnetic (AFM) quantum critical point 
(QCP), following the relation 1/Ti oc T/(T - T N ) 1/2 with T N = 0.06 K in the range of T = 
1.3 — 25 K. The onset of either the spin-density-wave (SDW) or charge-density-wave (CDW) 
order at To = 0.9 K, that is, of the first order, is evidenced by a broadening of the NQR 
spectrum and a marked reduction in 1/Ti just below To. The /-electron-derived correlated 
band realized in CeOs4Sbi2 is demonstrated to give rise to the novel phase transition on the 
verge of AFM QCP. 
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Filled-skutterudite compounds RcT4Pni2 (Rc = rare 
earth; T = Fc, Ru or Os; Pn = P, As or Sb) have at- 
tracted much attention because of their rich physical phe- 
nomena that are not yet fully understood. For instance, 
PrFe4Pi2 shows a heavy fermion (HF)-like behavior with 
a large mass of m* ~ 70 m e under a magnetic field (H) 1 
and undergoes an anomalous transition at a temperature 
T = 6.4KatiT = 0, indicative of a quadrupolar order- 
ing. 2-4 Note that the large C/T value for this compound 
is due not only to low-energy degrees of freedom such 
as either magnetic or quadrupolar fluctuations, but also 
to the Schottky anomaly originating from some low-lying 
CEF splitting. PrOs 4 Sbi2 is the first Pr-based HF super- 
conductor to reveal the large jump in the specific heat at 
T c = 1.85 K, the slope of the upper critical field near T C) 
and the electronic specific heat coefficient 7 ~ 350 — 500 
mJ/K 2 mol in the normal state. 5 The CEF-energy scheme 
of PrOs4Sbi2 has recently been determined to be the sin- 
glet Ti ground state accompanying the very low lying ex- 
cited state T^, 6 ' 7 and it was suggested that the HF-like 
behavior exhibited by PrOs4Sbi2 may be relevant to a 
very close energy level between Ti and T 4 . Thus, in the 
Pr-based compound, the 4/ 2 -derived CEF effect plays an 
important role in their rich physical phenomena. 

In contrast, most Ce-based filled-skutterudite com- 
pounds show semiconducting behavior, on the basis 
of which they are called hybridization-gap semiconduc- 
tors. A series of CeT4Pi 2 (T = Fe, Ru or Os) com- 
pounds have a hybridization gap of 400~1500 K, and 
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as the lattice constant increases, the energy gap de- 
creases. 8 ' 9 CeT4Sbi2 compounds show semimetallic be- 
havior. CeOs4Sbi2 is, on the one hand, suggested to 
exhibit Kondo insulating behavior with a large specific 
heat coefficient, 7 ~ 92 mJ/K 2 mol, and a very small 
gap of about A/fee ~ 10 K at the Fermi level. 10 Filled- 
skutterudite compounds form a unique cubic structure 
(space group: im3, T^ No. 204). The rare-earth ion in 
this structure is surrounded by twelve Sb atoms that are 
strongly hybridized with nearly localized 4/ electrons. 
Band calculations on CeOs4Sbi2 revealed a unique band 
structure where no band gap is formed at the Fermi 
level. 11 Electrical-resistivity measurement under pres- 
sure suggests that a hopping conductivity mechanism is 
realized at low temperatures. 12 Furthermore, CeOs4Sbi2 
was reported to exhibit an anomaly at approximately 
T = 1 K. The entropy release below this temperature 
is not sufficiently large to be considered to arise from 
some intrinsic phase transition, but is indicative of some 
impurity effect. 10 On the other hand, measurements of 
specific heat, magnetoresistance and the Hall effect un- 
der H revealed that a phase transition takes place. 13,14 

Here, we report on a novel phase transition in 
CeOs4Sbi2 that emerges at T = 0.9 K on the verge of an 
antiferromagnetic (AFM) quantum critical point (QCP) 
via measurements of the nuclear-spin lattice relaxation 
rate 1/Ti and nuclear quadrupole resonance (NQR) spec- 
trum of Sb nuclei. 

Single crystals of CeOs4Sbi2 were grown by the Sb- 
flux method. 1 The Sb-NQR measurement on the single 
crystal, which was crushed into powder, was performed 
using the conventional spin-echo method at H = and 
in the range of T = 0.2 — 300 K using a 3 He- 4 He-dilution 
refrigerator. 

For CeOs4Sbi2, five NQR transitions at 4.2 K are 
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Fig. 1. 121 Sb- and 123 Sb-NQR spectra for CeOs 4 Sbi2 at T = 4.2 
K. 



shown in Fig.l. These five NQR transitions come from 
the two Sb isotopes 121 Sb and 123 Sb with natural abun- 
dances 57.3% and 42.7% and nuclear spins I = 5/2 and 
7/2, respectively, and with a nuclear quadrupole mo- 
ment ratio of 123 Q/ 121 Q ~ 1.361. Thus, two and three 
NQR transitions are observed for 121 Sb and 123 Sb, re- 
spectively. 15 The NQR Hamiltonian is described as 



H Q = h -fWl-I{I 



+ !-)}, (1) 



where I is the nuclear spin, vq = 3e 2 qQ /2I(2I — l)h 
is the nuclear quadrupole frequency, and i] is the asym- 
metry parameter defined as r\ = (V xx — V yy )/V zz . Here, 
V xx , V yy , and V zz are the components of the electric field 
gradient (EFG) tensor. The values of respective 121 Sb 



and 123 

123 



Sb NQR frequencies, 121 i/ Q ~ 43.861 MHz and 



vq ~ 26.630 MHz, are estimated together with 77 ~ 
0.463. These two NQR frequencies are due to the differ- 
ence in nuclear quadrupole moment Q, i.e., 123 Q/ 121 Q ~ 
1.361. CeOs4Sbi2 and PrOs4Sbi2 have the same crystal 
structure. It is then expected that the values of vq and 
rj are almost equivalent to those in PrOs4Sbi2 reported 
previously. 16 The full width at half maximum (FWHM) 
A/ in the NQR spectrum is as small as ~ 65 kHz at 4.2 
K, indicating the high quality of the sample. 

T x is measured for both the 123 Sb-2i/ Q (±3/2 <-> ±5/2) 
and 121 Sb-2i/Q transitions above 1.4 K. For NQR-Ti, note 
that either a magnetic or a quadrupole relaxation process 
or both are possible as its origin. In case of the magnetic 
process, 1/Ti should be proportional to the square of 
the nuclear gyromagnetic ratio (7^) that is scaled to the 
nuclear magnetic moment. This is the present case be- 
cause the values of (1/Ti)/7n 2 are equivalent for 121 Sb 
and 123 Sb, where 121 7N = 10.2549 and 123 7N = 5.5530 
MHz/T, respectively. 15 Ti was uniquely determined by 
a theoretical curve for the recovery of nuclear magne- 
tization m(t) in which the asymmetry parameters were 
incorporated. 17 In the case of I = 7/2 and 77 = 0.46, the 
recovery curve is given by 
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Fig. 2. The temperature (T) dependence of 1/Ti for 123 Sb-2i^Q 
transition at H = 0. The open (solid) triangles show l/Ti in 
the T decrease (increase) process. 1/Ti does not show the T 
hysteresis above 1.4 K (open circles). The solid line is the T 
dependence of 1/Ti calculated by assuming a rectangular shape 
of density of states (DOS) with a bandwidth 2D = 3000 K, a gap 
size of 2A = 320 K and an amount of residual DOS Nq/N tcb = 
0.3 (see text) as shown in the inset. 



+0.9065 exp 



-17.207* 
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The T dependence of 1/Ti for 123 Sb-2z/Q is indicated in 
Fig. 2. The 1/Ti in the range of T = 200 - 300 K de- 
creases exponentially, and then 1/Ti is well fitted by an 
activation type of formula 1/Ti ~ exp(— A//cbT) with 
a gap value of A/fes ~ 330 K. The gap value coincides 
with the energy separation A/feg ~ 327 K between the 
CEF ground state Ti and the excited state T$. Here, 
A/ks ~ 327 K is cited from a measurement of mag- 
netic susceptibility. 10 Such 1/Ti is, however, observed in 
the related compound CeRu4Sbi2, where the magnetic 
susceptibility measurement does not reveal a clear CEF 
effect. 18 ' 19 Also, recent ultrasound measurement has not 
detected a clear CEF effect on CeOs4Sbi 2 . 20 Therefore, 
this exponential decrease in 1/Ti may not be related to 
the CEF effect. In heavy-fermion systems, 1/Ti is domi- 
nated by /-electron-derived spin fluctuations through the 
hybridization between / and conduction electrons. 
In general, 1/Ti is given by 21 



1 

Ti 



27nfeeT 
(7cft) 2 



w 



(3) 



where A q is the q-dependent hyperfine coupling constant, 
and x± is the perpendicular component of the imaginary 
part of dynamical susceptibility. Provided that the q de- 
pendence of spin fluctuations is not so large, eq. (3) can 
be rewritten as 



1 

— oc T 
Ti 



N eff( E ) 



(4) 



where f{E) is the Fermi distribution function and 
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Fig. 4. The T dependence of 1/T\T at low temperatures. Open 
and closed triangles show the T dependence of 1/T\T upon cool- 
ing and heating, respectively. 



Fig. 3. The T dependence of 1/T\T where open circles corre- 
spond to the data of 1/Ti (open circles) in Fig. 2. The data 
in the range of T = 1.3 — 20 K are fitted by the relation 
1/Ti oc T/(T - Tn) 1 / 2 with T N = 0.06 K as shown by the 
solid curve. Note that this relation is consistent with the SCR 
theory for three-dimensional itinerant weakly antiferromagnetic 
metals 25-27 The inset shows a 123 (TiT) 2 vs T plot where the 
solid line corresponds to the relation (TiT) 2 oc (T— Tn), ensuring 
the theoretical prediction in terms of the SCR theory mentioned 
above. 

N e ff(E) is the rcnormalizcd quasiparticle density of 
state (DOS) at the Fermi level. For "Kondo" insula- 
tors reported to date, 1/Ti is described by a pseudogap 
model induced by the hybridization between / and con- 
duction electrons. In CeRhSb and CeNiSb, which form 
orthorhombic structures, 1/Ti is explained by the V- 
shaped gap model. 22 ' 23 In Ce3Bi4Pt 3 , which forms a cu- 
bic structure, 1/Ti that decreases exponentially is ex- 
plained by assuming an isotropic gap model with a rect- 
angular DOS. 24 The difference in the gap structure for 
the "Kondo" insulators may be related to their crystal 
structure. In this context, for CeOs4Sbi2, which forms 
a cubic structure, the isotropic gap model with a rect- 
angular DOS is expected to be applied for interpreting 
1/Ti. The solid line in Fig. 2 shows 1/Ti calculated by 
the model DOS illustrated in the inset. Here, the fitting 
parameters are 2D = 3000 K and 2A = 320 K with a 
finite residual DOS No/N les = 0.3 at the Fermi level. 
This model explains the 1/Ti above T ~ 90 K. If the 
residual DOS at the Fermi level existing inside the pseu- 
dogap is responsible for the relaxation process at a low 
T, a behavior in which TiT = constant would be ex- 
pected at temperatures lower than A/ks- Actually, al- 
though a behavior in which T\T = constant is observed, 
such behavior is only valid in the range T = 25 — 50 K. 
Unexpectedly, \jT\T increases upon cooling below 25 
K as shown in Fig. 3. This enhancement in 1/T±T has 
not been observed for other "Kondo" insulators. It is re- 
markable that 1/Ti behaves as 1/Ti oc \/T below T = 25 
K. The self-consistent renormalization (SCR) theory for 
spin fluctuations in itinerant AFM metals has been ap- 



plied to the case of HF systems near an AFM QCP. 25 27 
In this case, it was predicted that 1/Ti obeys the relation 

y^ t 4^-^jt=% (5) 

because of the Currie- Weiss law of staggered susceptibil- 
ity Xq(T) oc 1/(T— Tn) with an AFM ordering temper- 
ature Tn- Therefore, note that 1/Ti oc y/T is valid at 
T N ~ in the vicinity of an AFM QCP. In order to fur- 
ther inspect this unique T dependence, (TiT) 2 is plotted 
as a function of T in the inset of Fig. 3. From this plot, 
it is ensured that the relation (TiT) 2 oc (T — Tn) is valid 
with Tn = 0.06 K in the range T = 1.3 — 25 K. It is 
interesting that CeOs4Sbi2 is near the AFM QCP and 
an estimation of "Tn = 0.06 K" in Xq(T) is indicative 
of a possible onset of some AFM order at low temper- 
atures. As a matter of fact, suggesting the appearance 
of some anomaly, the recovery curve below 1.4 K is not 
described by eq. (1). Therefore, a short component in 
Ti is plotted below 1.4 K in Fig. 4. With further de- 
creasing T, a marked decrease in 1/Ti is observed below 
T = 0.9 K. Figure 4 shows the detailed T dependence 
of 1/TiT around T = 0.9 K. Unexpectedly, this T x re- 
sult reveals the onset of some phase transition at T=0.9 
K, which is higher than "Tn = 0.06 K" estimated from 
Xq(T). The sudden decrease in 1/TiT is due to the open- 
ing of a gap at the Fermi surface, caused by the onset 
of phase transition. The data are consistent with such 
an activation type of behavior as 1/Ti oc exp(— A/k^T) 
with a gap size of A/fee ~ 1-83 K. It should be noted 
that 1/TiT undergoes clear hysteresis upon cooling and 
heating. The hysteresis in 1/TiT suggests a first-order 
transition emerging at T = 0.9 K, which is also corrob- 
orated by the measurement of the T dependence of the 
NQR spectral shape as shown later. 

Figure 5(a) indicates the T dependence of the NQR 
spectrum at the 123 Sb-2z/Q transition below 4.2 K. Above 
0.9 K, it exhibits a Lorentzian shape with a small value 
of Af = 65 kHz. With decreasing T below 0.9 K, it starts 
to broaden as clearly seen in the figure, and at the same 
time, 1/Ti is markedly reduced. Note that the tail in the 
NQR spectrum is significantly increased and makes the 
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Fig. 5. (a) The T dependence of the NQR spectrum of 123 Sb- 
2uq transitions, (b) The T dependence of FWHM upon cooling 
(closed circles) and heating (open squares). 



NQR spectral shape deviate from its Lorentzian form. 
This result may be ascribed to either the appearance of 
an internal field that makes its magnitude distributed at 
the Sb site or a distribution of i>q. Then, such a charac- 
teristic spectral shape below T = 0.9 K points to either a 
wide distribution of the internal field at Sb if a transition 
were a spin density wave (SDW), or a distribution in the 
electric field gradient associated with a charge density 
wave (CDW), both presumably triggered by the nesting 
effect of the Fermi surface. Furthermore, the T depen- 
dence of FWHM in the NQR spectrum shows hysteresis 
upon cooling and heating as shown in Fig. 5(b). This re- 
sult also suggests that the first-order transition begins to 
take place below T = 0.9 K, consistent with the hystere- 
sis observed in the T dependence in 1/Tj. The present 
results on CeOs4Sbi2 have revealed that the phase tran- 
sition near the AFM QCP is of the first order. This find- 
ing deserves a theoretical study on the physics behind 
the AFM QCP in /-electron-derived correlated bands 
realized in HF systems in general. 

In summary, we have shown from the measurements of 
1/Ti and the NQR spectrum of Sb nuclei that CeOs4Sbi 2 
is a unique Kondo semiconductor. CeOs4Sbi2 has a hy- 
bridization gap similar to other Kondo semiconductors, 
however, there exists a large residual DOS within the 
gap. The enhancement of spin fluctuations below 20 K 
indicates that CeOs4Sbi2 is closely located to the AFM 
QCP. Even in such a case, the remarkable finding is that 
either a SDW or CDW order, which is of the first order, 
takes place at T = 0.9 K, possibly triggered by the nest- 
ing effect of the Fermi surface. The present result, we 
believe, deserves a theoretical study on the physics be- 
hind the AFM QCP realized in HF systems in general. 

note added in proof. After the submission of this pa- 
per, we became aware of an important theoretical work 
by T. Ohashi, A. Koga, S. Suga, and N. Kawakami (Phys. 
Rev. B70 (2004) 245104). This work has suggested that 
if some state were induced inside the Kondo gap, the an- 
tiferromagnetic correlation might be developed, leading 
to an SDW instability at low temperatures as observed 
in the present study. 
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